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Abstract 
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systems related to generalized statistics. We discuss several fermionic, 
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1 Introduction 



The Bethe ansatz is a technique which is based upon the quantum mechanical 
description of a many particle system by a wave function. The interaction of the 
individual particles is assumed to be relativistic, short-range and characterized by 
a factorizable scattering matrix. The boundary condition for the many particle 
wave function leads to what is commonly referred to as the Bethe ansatz equation, 
which provides the quantization condition for possible momenta of this system. 
Taking the thermodynamic limit of this equation, that is taking the size of the 
quantizing system to infinity, leads to the so-called thermodynamic Bethe ansatz 
(TBA). The origins of this analysis trace back to the seminal papers by Yang and 
Yang m and the technique has been refined and applied to different situations 
in numerous works [2-6] thereafter. The TBA constitutes an interface between 
massive integrable models and conformal field theories. One may extract different 
types of information from it, where the ultraviolet behaviour, i.e. ultimately the 
central charge of the conformal field theory, is the most accessible. 

In the derivation of the TBA-equation the underlying statistical interac- 
tion (also called exclusion statistics) is usually taken to be either of bosonic 
or fermionic type. Seven years ago Haldane [|7| proposed a generalized statistics 
based upon a generalization of Pauli's exclusion principle. This type of statis- 
tics (anyonic) has many important applications, in particular in the description 
of the fractional quantum Hall effect M. The main purpose of this manuscript 
is to implement systematically the Haldane statistics into the analysis of the 
thermodynamic Bethe ansatz. Hitherto, attempts in this direction || 10 



were 



mainly based on the consideration of particular statistical interaction which only 
involves one species, like the Calogero-Sutherland model . Our approach will 



cover a general choice of statistical interaction described by some, in general non- 
diagonal, matrix g^. We put particular emphasis on the ultraviolet region, which 
corresponds to the high temperature regime. We also continue our investigation 
started in |Ll| and clarify the role of the anyonic S-matrix (p5| ) in the context of 



the TBA. 

We formulate a macroscopical equivalence principle in the sense that the 
macroscopical nature of a multi-particle system is only governed by the com- 
bination flT3| ) of the dynamical- and statistical interactions. This means in par- 
ticular that two multi-particle systems differing on the microscopical level, i.e. in 
the S-matrix, may be made macroscopically equivalent by tuning the statistical 
interaction. 

Our manuscript is organized as follows: In section 2 we recall the derivation of 
several thermodynamic quantities from Haldane statistics. In section 3 we derive 
the thermodynamic Bethe ansatz equation for a multi-particle system in which 
the statistical interaction is governed by Haldane statistics. In section 4 we ar- 
gue that certain multi-particle systems may be transformed into macroscopically 
equivalent systems by tuning the statistical- and the dynamical interaction. In 
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section 5 we demonstrate that certain scattering matrices, leading to equivalent 
thermodynamical systems differ only by CDD-ambiguities and comment on the 
ambiguity in choosing a particular statistical interaction. In section 6 we discuss 
the ultraviolet limit of the generalized TBA-equation. In section 7 we derive 
Y-systems related to generalized statistics. In section 8 we illustrate our general 
statements by some explicit examples. Our conclusions are stated in section 9. 



2 Thermodynamics from Haldane Statistics 

The object of our consideration is a multi-particle system containing I different 
species confined to a finite region of size L. We denote by rii the number of 
particles, by iVj the dimension of the Fock-space related to the species i and by 
di the number of available states (holes) before the nj-th particle has been added 
to the system. When treating bosons the number of available states naturally 
equals the total dimension of the Fock space, i.e. iVj = <ij, whereas when treating 
fermions there will be restrictions due to Pauli's exclusion principle and one has 
Ni = di + rii — 1 . With these relations in mind the total dimension of the Hilbert 
space may be written for both cases as 

w Md i + ni -l)\ 

t\ mKdi-iy. { ' 

Conventionally one employs in (|l|) the Fock-space dimension rather than the 
number of available states before the n^-th particle has been added to the system. 
However, besides being a unified formulation, equation ([!]) has the virtue that it 
allows for a generalization to Haldane statistics. 

By introducing a statistical interaction g^, Haldane @ proposed the following 
generalized Pauli exclusion principle 

^ = - g - ■ (2) 

Relation (|2]) means that the number of available states should be regarded as a 
function of the particles present inside the system. In this proposal a statistical 
interaction between different particle species is conceivable. In the bosonic case 
there will be no restriction such that g^ = 0, whereas in the fermionic case the 
number of available states reduces by one if a particle is added to the system, 
hence one chooses gij = 5ij. One assumes |7], ||, that the total dimension of the 
Hilbert space is still given by (|I|), where the quantities involved are related to 
each other by @. 

We now want to analyze the multi-particle system in its thermodynamic limit, 
that is we let the size of the confining region approach infinity, L — > oo. It is then 
a common assumption that the ratio of the particle (hole or state) numbers over 
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the system size L remains finite. For the fraction of particles (holes or states) 
of species % with rapidities between 9 — A6/2 and 9 + A8/2 it is convenient to 
introduce densities 

ANi = Pi {6)A6L (3) 
An t = tf(0)A0L (4) 
Ad, = p?(0)A0L . (5) 

The rapidity 9 parameterizes as usual the two-momentum p = m (cosh#, sinh#). 
Integration of the generalized Pauli exclusion principle (|2]) then yields a relation 
between the different types of densities 

/hV)=AV) + £,9*4(0) ■ (6) 

The constant of integration has been identified with the Fock space dimension. 
The reason for this identification is based on fact that in this way one recovers 
for gij = and <?jj = 5ij the usual bosonic and fermionic relations, respectively. 
Notice that in the finite case the constant is in general slightly different from 
TV,, for instance for fermionic statistics it has to be chosen as iVj + 1 in order to 
recover the relation Ni = d% + n, — 1. However, in the thermodynamic limit this 
difference is negligible. 

Let us remark, that equation @ suggests from a physical point of view that 
gij should be non-negative. However, if there are additional symmetries then this 
requirement can be weakened. For instance for two conjugate particles, say % and 
% (see section 8 for examples), one naturally assumes that p\{&) = pl(0), such 
that only the combination + has to be taken non-negative. 

We are now in the position to construct further thermodynamic quantities. 
First of all we may sum up all contributions from occupied states in order to 
obtain the total energy 

i 00 

E[p r ] = Lj2 J d9p r i{Q)™-i cosh 9 ■ ( 7 ) 

i=1 -oo 

Furthermore, we obtain from ([]]), upon using Stirling's formula Inn! ~ nlnn, 
(H)-(H) an d ©, the entropy S = k\nW as a functional of the particle- and 
Fock-space density 



i 

S[p,p r ]=kLY: J d6 



Pi + h ijP ^\ ( Pi + h ijP y 



1 \Pi-9ijPjJ V Pi 



We introduced here = 5^ — g^ and use the sum convention (sum over j) to 
avoid bulky expressions, k is Boltzmann's constant. According to the fundamen- 
tal postulates of thermodynamics the equilibrium state of a system is found by 



3 



minimizing the free energy F. Hence, keeping the temperature constant we may 
obtain the equilibrium condition by minimizing F [p, p r ] = E [p r } — TS [p, p r ] 
with respect to p r . The equilibrium condition reads 

6F 5E ^5S m ' 8S 5 Pj 

h\ $pi dpi px^pjSpl 

So far we did not provide any information about the admissible momenta in the 
system, which are restricted by the boundary conditions. 



3 Thermal Equilibrium with Boundaries 

Boundary conditions may be accounted for by the Bethe ansatz equations. Recall 
that the Bethe ansatz equation is simply the equation which results from taking 
a particle in the multi-particle wave function on a trip through the whole system 
0. The particle will scatter with all other particles in the system, described by 
a factorizable S-matrix, such that 

i 

exp(iLrriiSmh.9i) TT Sij(9i — 9j) = 1 , (10) 

has to hold for consistency |l|. To simplify notations we may assume here that the 
scattering matrix is diagonal, such that the subscripts only label particle species. 
This set of transcendental equations determines which rapidities are admissible 
in the system due to the quantization as a result of restricting the size of the 



system. Taking the logarithmic derivative of the Bethe ansatz equation (10) and 
employing densities as in (||) and (§) one obtains (see for instance ||, or 
for more details) 



^7n i coahe + J E(<Pij*Pj)(0)=Pi(0) ■ (11) 



2tt 

Here we introduced as usual the notation <fiij(9) = —i4n In 5^(0) and denote the 
convolution by (/ * g) (6) := 1/(2-71") / d9' f(9 — 6')g{9'). The equilibrium condition 
(|9D together with (|TT|) yields the desired thermodynamic Bethe ansatz equation 
for a system in which statistical interaction is governed by Haldane statistics 

1 1 
— m l cosh^ = ln(l + x l (^)) + E(^'* ln (l + a; 7 1 )) ( d ) ■ ( i2 ) 

We use here the abbreviations Xi (9) := p\ (9) j p\ (9) and 

$ ij (e):=(p ij (e)-2irg ij 5(e) . (13) 
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In the derivation of ([12]) we assumed that = gji. We assume that Xi (9) is 
symmetricQ in the rapidity throughout the manuscript. In general one is only 
able to solve equations (|T2| ) numerically as we demonstrate in section 8. In the 
formulation of the TBA-equation of bosonic or fermionic type, it is common to 
introduce here as an additional quantity the so-called pseudo-energies Ei{9). In 
general one may employ 

ln(l + x i (9))-^2g ij \n(l+xj 1 (9))^e i (9) , (14) 
i=i 

for this purpose. For the finite case the same relation was obtained by Wu 0, 
there however, the quantity X{(9) has a slightly different meaning. Clearly it is not 
possible to give a general solution of fli~4D . However, for the bosonic and fermionic 
case it is solved easily, we obtain X{(9) = exp(£j(0)) — 1 and %i{9) = exp(£j(0)). 
Substitution of these solutions into equations fljjD turns them into the well-known 
TB A-equations of bosonic and fermionic type , respectively. With the ratios 
of the related particle- and Fock-space densities 

Pl {9) exp(£^))Tl 1 ) 

become the usual Bose-Einstein (upper sign) and Fermi-Dirac (lower sign) distri- 
butions. From our point of view it does not seems to be necessary to introduce 
pseudo-energies and one should rather view Xi (9) = p\ (9) / p[ {9) as a more fun- 
damental entity. In addition one avoids the problem of solving (|14T). From a 
numerical point of view however, it appears sometimes useful to formulate the 
TBA-equation in different variables. 

We now substitute back the equilibrium condition into the expression for the 
free energy and obtain together with the general expressions for the total energy 
(0) and the entropy (||]), the generalized thermodynamic Bethe ansatz equation 
©and© 

LkT 1 7 

F(T) = / d9m i cosh 9 ln { 1 + x 7 1 ( )) ■ ( 16 ) 

The relation between the free energy and the finite size scaling function is well- 
known to be c(T) = —6F(T)/(7rLT 2 ) Jl4|]. As usual we now identify the tem- 
perature with the inverse of one radial size of a torusi T = 1/r and choose now 
Boltzmann's constant to be one. Then 

c ( r ) = % Y, m *] d9cosh9ln ( 1 + x ^ 1 ( 9 )) ■ ( 17 ) 

n i=1 o 

* Taking into account that $ij (9) is symmetric in 8 (due to (|l])), this is equivalent to the 
assumption that jl^ ) may be solved iteratively (see e.g. |fl3|). 

tFor more details on the physical picture see for instance g [| § [lj. 
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Once more for Xi(9) = exp(£j(0)) — 1 and Xi(8) = exp (£;(#)) we recover the well- 
known expressions for the scaling functions of bosonic and fermionic statistics. 
In the ultraviolet limit the scaling function becomes the effective central charge 
of a conformal field theory |14[], i.e. lim r ^o c(r) = c e // = c — 24h'. Here c is the 
usual conformal anomaly and h! denotes the lowest conformal dimension ]T3 



To summarize: For a given statistical interaction (0) and dynamical interac- 
tion described by a factorizable scattering matrix, we may solve in principle the 
generalized thermodynamic Bethe ansatz equation ([T2|) for Xi(6). This solution 
together with the knowledge of the mass spectrum of the theory, allows (up to 
a one dimensional integral which may always be carried out by simple numerics) 
the calculation of the entire scaling function. 



4 Equivalent Multi-Particle Systems 

From a thermodynamic point of view, systems which have the same expressions 
for the free energy show the same behaviour. Hence multi-particle systems which 
possess the same scaling function over the entire range of the scaling parameter r 
are to be considered as equivalent. This implies that two multi-particle systems 
are equivalent, if the mass spectra and the quantities 3>y are identical. As a 
consequence of this we may achieve that two systems equal each other from a 
macroscopical point of view, despite the fact that they involve different scattering 
matrices describing the dynamical particle interaction. The apparent difference 
can be compensated by a different choice of the underlying statistical interaction. 
Considering for instance the well-known TBA-equations of bosonic- and fermionic 
type (see for instance J2|) involving (p\j{9) and <p{j{9), respectively, we observe that 

dm=vU0) + ^5 t3 5{6) (18) 

transforms both equations into each other. 

The question of whether relations between two scattering matrices leading to 
(|T8|) are at all conceivable immediately comes to mind. Such relations emerge in 
several places. For example in || the authors assume 

S[ 3 {9) = Sij(e) exp(-27ri<f y e(0)) . (19) 

The function ©(#) was taken to be the usual step-function, with the property that 
O(0) = 1/2, such that Sij(9) and S'^9) only differ at the origin of the complex 
rapidity plane. This modification was necessary in order to extend the validity of 
certain identities ((2.18) in || leading to eq. ( fTO|) below) involving the scattering 
matrices of ADE-afline Toda field theories (see for instance ||16||), i.e. Sij(9), also 
to the origin. In particular, in some special cases these identities become equal 
to the bootstrap equations. A further motivation to introduce (|T9"D was to derive 
so-called Y-systems proposed by Zamolodchikov (last ref. in 0) for ADE-afline 
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Toda field theories. In this case the system which involves Sij(9) with fermionic 
statistical interaction is equivalent to the system involving S'^iO) with bosonic 
statistical interaction. 

As a further example we may also consider a system in which the dynamical 
interaction is described by the scattering matrix S{6) = — exp(— iir\e(9)), where 
A denotes the coupling constant and e{9) = 0(0) — 0(— 9), of the Calogero- 
Sutherland model with a statistical interaction of bosonic type, i.e. g = 0. 
This system is equivalent to a system with a constant S-matrix, e.g. S'(9) = — 1, 
and a statistical interaction of the form g' = A. 

In general we have the following equivalence principle. Multi-particle systems, 
involving quantities such that the relation 

tp' ij (9) = t Pij (9)+27r(g' ij -g ij )S(9) (20) 

holds and the masses for the same species are identical, are thermodynamically 
equivalent. 



5 Microscopical Ambiguities 
5.1 On CDD- Ambiguities 

We now want to elaborate on the question in which sense the scattering matrices 
leading to (|20|) differ. The analysis of analytic properties of the scattering matrix 
leads to a set of consistency equations, which have to be satisfied by any S- 
matrix related to integrable models in 1+1 dimensions. These equations are so 
restrictive, that they determine the S-matrices of a particular model up to what 



is usually referred to as CDD-ambiguities [17], [18| . In this section we argue that 
certain S-matrices related to each other in such a way, that they lead to (20), 
differ precisely by such ambiguities. 

A scattering matrix Sij(9) is usually assumed to be a meromorphic function 
in the strip < lm.9 < it,— oo < lie 9 < +oo. This region is considered as 
physical, meaning that all singularities occurring in this sheet acquire a physical 
interpretation. The scattering matrix in fllQ|) is usually regarded as the one which 
results from an analysis of the so-called bootstrap equations 

SiMSiA-o) = i (2i) 

S l3 (9)S- l3 (9-m) = 1 (22) 
Suie + i^Siiie + iri') = S fk (9) . (23) 



Equation (^T|) is a result of unitarity and analytic continuation, (22) a conse- 
quence of crossing invariance and (^) (strictly speaking this is the boostrap 
equation) expresses the factorization property for the fusing process i + j — > k. 
The so-called fusing angles rj, r( are specific to each model depending on the mass 
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spectrum. It was found in [18| that the most general solution to these consistency 



equations will always be of the form 

tanh \{9 + a) 
i^tanh^-a) ' 1 J 

where A is a set of complex numbers which characterizes a particular model. 
When particles are not self-conjugate one should replace tanh by sinh in 



There is however the freedom to multiply these expressions with a so-called CDD- 
ambiguity [T7]] (also of the form (JMp, but related to a different set A 1 ), which 
satisfies by itself all the consistency requirements without introducing any addi- 
tional poles into the physical sheet. A well-known example for such an ambiguity 
are for instance the coupling constant dependent blocks. One seeks solutions of 
the general form fl2~4|) involving two sets A and A' such that the product related 
to the set A (the so-called minimal S-matrix) already closes the bootstrap and 
accounts for the whole particle spectrum independent of the coupling constant. 
Then the additional factors related to A' constitute a CDD-ambiguity depending 
on the coupling constant f3 in such a way, that in the limit (3 — > (in theories 
which admit duality also (3 — > oo) the S-matrix becomes free, that is one. 

As a particular case of two S-matrices related to each other such that they 
may lead to an equation of the type (p0|) we will now consider 



= Sij(8) exp(-27rzA+.9(Re#) - 2nA^0(- Re 9)) , (25) 

which we discussed in [jP2| . Similarly as in || we choose 0(0) = 1/2. Here the 
A^ are related to the asymptotic phases of the S-matrix Sij(9) 



lim £y(0) = exp(27riAg . (26) 

He y^±oo 



The asymptotic phases are well defined, since the limit Re 9 — > ±oo of ( P^ ) does 
not depend on the imaginary part of 9. This property is in particular needed 
to obtain (P9"|). The transformation fl2"5j) compensates the asymptotic phases and 



creates a non-trivial phase at 9 — > 0, i.e. the anyonic situation. It was argued (19 
that from a physical point of view, S-matrices which possess a non-trivial asymp- 
totic phase should be regarded rather as auxiliary objects. Scattering matrices of 
the type ( ^5|) should be considered as the genuine physical quantities, since they 
lead to the correct physical properties, including the exchange statistics. 



We shall demonstrate that the system of equations fl21|)-(|2~3D does not have 
to be altered for the anyonic matrix S'^{9). For this purpose we will first derive 



some properties concerning the phases. Taking the limit Re 9 — > oo in (pl|) 
yields immediately several relations between the asymptotic phases 



A+ + A7. = n (27) 
AJ + AJ = n' ± (28) 

v ij - ^ik 



A± + A± = Aj + < , (29) 
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with n,n±,n± being some integers. Using now these relations for the phases 
together with the fact that Sij{6) satisfies the consistency equations (|21p-(|23D, 
it is straightforward to derive the related equations for the anyonic scattering 
matrix S'^iO) 

s' l0 {e)s' l0 {-e) = 1 (so) 

S'^S'^e-m) = 1 (31) 
S^e + i^S'^e + irj') = S' fk (9) . (32) 

Comparing and 5 , - J (6 ) ), we conclude, taking ( p7f) into account, that the 

additional factor in (p5|) has altered the behaviour at the imaginary axis of the 
rapidity plane only up to a sign. Thus Sij{9) and S'^iQ) differ by a CDD-factor. 

It is straightforward to generalize the previous argument to the case when the 
scattering matrix is non-diagonal, such that in addition to (pl]) - (p3|) one also has 
to satisfy the Yang-Baxter equation as a consequence of factorization. 

5.2 On statistical Ambiguities 

In this subsection we recall the argument which leads to a particular choice of the 
statistical interaction. In general one considers the value of the scattering matrix 
at 9 = in order to deduce the symmetry properties of the Bethe wave function. 
Then together with the a priori (or e.g. from a Lagrangian) knowledge of the 
nature of the particles one deduces the statistical interaction. For example for 
Su{0) = — 1 the Bethe wave function is antisymmetric. If in this case one describes 
bosons, one is forced to choose a fermionic statistical interaction. However, this 
way of arguing seems somewhat ambiguous as a simple example demonstrates. 
Considering for instance a system in which the dynamical interaction is described 
by the affine Toda scattering matrix, one usually selects @] fermionic statistics 
due to the fact that 5*^(0) = (— l) 5ij . In the limit (3 — > (because of strong- 
weak duality one may also take (3 — *> oo) the theory becomes free and one obtains 
Sij(8) = 1. Using now the same arguments one has to deduce, from the symmetry 
of the Bethe wave function and the fact that one still describes bosons, that the 
statistical interaction has to be bosonic. Concerning the scaling function there 
is no problem here since in both cases we obtain the same ultraviolet limit. 
However, there is a change in the statistics. The apparent paradox is resolved 
by making use of the macroscopical equivalence. Using the scattering matrix 
(|19|) instead one has now SV, (0) = 1. Making then use of the bosonic-fermionic 
transformation ( Jl8| ) we obtain by means of the same arguments a unique, that is 
bosonic, statistical interaction for the entire range of the coupling constant (3. 

Scattering matrices related to each other as in (|2"o] ) play a distinguished role in 
this context since they only differ by a CDD-ambiguity as argued in the previous 
section. Considering such expressions the relation between the corresponding 
quantities ^{9) and <fij{9) are fixed. In this case we achieve thermodynamical 
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equivalence by demanding 



g' ij -g ij = Ar j -At j . (33) 

We may consider a few examples. For instance in the case of affine Toda field the- 
ory we have A^ = ±5^/2 0, such that gL — gij = Sij, which is a transformation 
from fermionic exclusion statistics to bosonic one. In case we only consider the 
minimal part of the scattering matrix for ADE-affine Toda field theory we obtain 
from the asymptotic behaviour observed in f§], i.e. A^ = ±(5^/2 — (C -1 )^), 
that is g[j — = C~ l I , where C denotes the Cartan matrix and I the incidence 
matrix of the related Lie algebra. Assuming now that g^ is fixed by the argu- 
ments presented above, we obtain for the later case an interesting expression for 
the statistical interaction g'^ in terms of Lie algebraic quantities. Therefore one 
may formulate the generalized Pauli principle in this context in a Lie algebraic 
form. 



6 The ultraviolet Limit 

One of the interesting quantities which may be extracted from the thermodynamic 
Bethe ansatz is the effective central charge of the conformal field theory when 



taking the ultraviolet limit [14fl, i.e. r — > in fll7|) . The integral equation (|T 



simplifies in this case to a set of constant coupled non-linear equations 

i 

Hl + x i )=J2(N, lJ +g ll )\n(l + xj 1 ) , (34) 



3=1 



where = A^ — A^ . In this limit one may approximate rm $ cosh 9 in (|I2"D and 
(|1^) by exp6 rrrii/2 i. Taking thereafter the derivative of ( 12 ) we obtain, upon 



the substitution of the result into (|17|) , for the effective central charge 

-//^X>(t^-) • (35) 

Here L(x) = -\ J X dt + g] = £^ =1 g + \ lnxln(l - x) denotes Rogers 

dilogarithm In these definitions it is assumed that x takes its values between 
and 1, which in turn implies that all in ( p5[) are to be non-negative. This 
is in agreement with the physical interpretation of the Xi as ratios of densities. 
Once again with Xi = exp(si) — 1 and Xi = exp(si) we recover the well-known 
expressions for the bosonic and fermionic type of statistical interaction, respec- 
tively. 



*Of course these approximations rely upon certain assumptions (for details see for instance 
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Obviously the transformation properties discussed in section 4 also survive 
this limit process. As follows directly from ([34]) they read now 



= Nij + 9ij - ^ . (36) 

Apparently this condition is weaker than (|20|). It guarantees equivalence multi- 
particle systems only at the conformal point. The transformation from bosonic 



to fermionic statistics is compatible with equations (50) and (51) in |12| for the 
case b=l, where such transformations where obtained purely on the conformal 
level, that is by means of manipulations of certain characters. 

It is instructive to consider a few examples. For instance having the situation 
that the difference of the phases of the scattering matrix equals the negative of 
the statistical interaction, i.e. Nij = —gij we always obtain 

c e// = ^£L(l) = / . (37) 

An example for this situation is to consider a system involving the scattering 
matrix related to ADE-affine Toda field theory (recall that in this case Ny = —5ij 
0]) and choose the statistical interaction to be fermionic. 

An interesting structure emerges when considering a system in which the dy- 
namical interaction is described by a direct product of Calogero-Sutherland scat- 
tering matrices, i.e. Sij(6) = — exp(— «7rAi5ye(#)), and the statistical interaction 
is taken to be of fermionic type. Then ([34]) reduces to 



x} +Xi = {1 + Xi ) Xi . (38) 

The same equation is of course obtained when the statistical interaction is chosen 
to be bosonic and the coupling constants shifted by one. In some cases ( p8|) is 
solved easily analytically and we can employ fl3"5| ) to compute the effective central 
charge. We may decompose the effective central charge into contributions coming 
from different choices for the A«: c e // = Z)i=i c e//- For instance for Aj = 1 the 
solution of fl3"8p is Xi = (1 + v5)/2 and for Aj = —1/2 we obtain Xi = (\/E — l)/2. 
Then for Aj = —1/2 and A; = 1 we obtain, with the help of fl35|), 



= 5 and <" = ^(jT^H (39) 

respectively. As already mentioned the computation of the effective central charge 




by means of ( p4|) and (3~5]) is not always rigorously justified. However, the explicit 
analytic computations of the full scaling functions provided in the next section 
confirm these results. It is intriguing to note that we recover in this way the 



effective central charges for certain minimal models [ 15] . The case A = —1/2 



corresponds to the minimal model A4(3, 5), whilst A = 1 corresponds to the 
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Yang-Lee model, i.e. A4(2,5). Equation ([38]) may also be solved for different 
values of A and we may compute the effective central charge by means of fl35|) . 
The dependence of c e // on A for the one-particle contribution is depicted in 
figure 1. i This figure suggests immediately that one may find other minimal 
models as ultraviolet limits of ([T7|) for different values of A. However, apart from 
A = 0, —1, oo, the two values presented are the only possible choices for a model 
with solely one particle leading to accessi bleS relations for dilogarithms. This 
implies that other rational A do not lead to rational values of c e ff. 

7 Y-systems 

For some class of models it has turned out to be possible to carry out certain ma- 
nipulations on the TBA-equations such that the original integral TBA-equations 
acquire the form of a set of functional equations in new variables Yj (last ref. 
in [[|). These functional equations have the further virtue that unlike the orig- 
inal TBA-equations they do not involve the mass spectrum and are commonly 
referred to as F-systems. In these new variables certain periodicities in the ra- 
pidities are exhibited more clearly. These periodicities may then be utilized in 
order to express the quantity Y as a Fourier series, which in turn is useful to 
find solution of the TBA-equations and expand the scaling function as a power 
series in the scaling parameter. We will now demonstrate that similar equations 
may be derived for a multi-particle system in which the dynamical scattering is 
governed by the scattering matrix related to ADE-affine Toda field theories and 
the statistical interaction is of general type. 

We consider the modified version (in the sense of (|i~9|)) of the minimal part 
of the scattering matrix of ADE-affine Toda field theories. As was shown in || 
these S-matrices lead to the identity 

(17T \ / ITT \ r 

6 + ~h) + VlJ V~~h) = ^ Iik(fikj {e) ~ 2nIij6 {e) ' (40) 

where h denotes the Coxeter number, r the rank and I the incidence matrix of 
the Lie algebra. It is then straightforward to derive the U YZ— system" 

y '( e+ x) yi ( fl -f) = n^ (fl)/w - (41) 

§The plot is presented only for A > — 1 since eq. (|3^) possesses non- negative solutions only 
at this interval. 

^In general a dilogarithmic identity is called accessible if it is of the form yj i=1 L(yi) = kn 2 /6, 
with k being rational and yt algebraic H|. 
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The quantities Y (9) and Z (9) involve now the statistical interaction in the form 



Y t (6) = (l + x t (9))Y[(l + xj 1 (6))-^ and Z { (9) = Y, L (9) (1 + x~ l {6)) . 
i=i 

(42) 

Equations (fy]) follow upon first adding (|i~2l) at 9 + and 9 — ^ and subtracting 
I times (0) at 9 from the sum. Thereafter we employ the fact that the masses 
of affine Toda field theory are proportional to the Perron- Frobenius vector of the 

4sin 2 (7r/(2/i))m i 



Cartan matrix, i.e. Yfj=i Cij m j 



Then, with the help of 



(0) and flU), the equations (|T]) follow. 

In comparison with (|T^) the equations fl4"T| ) have already the virtue that they 
are simple functional equations and do not involve the mass spectrum. However, 
in order to solve them we still have to express Z in terms of Y or vice versa, 
which is not possible in general. However, once the statistical interaction is 
specified this may be achieved. For instance, for gr^ = Sij we obtain Zi = Yi + 1 
and recover the known fermionic F-system. In the bosonic case, i.e. gij = we 
obtain Z{ = Y- 2 / (Yi — 1) and for "semionic" statistical interaction with = Sij/2 



we obtain Zi 



1 + 4K 2 



8 Examples 

8.1 Ising model and Klein-Gordon Theory 

The most elementary examples which illustrate the features outlined above more 
concretely is simply to consider the Ising model (Ai-minimal affine Toda field 
theory) S(9) = —1 or the Klein-Gordon theory S(9) = 1. Then for both cases 
equation flT2| ) is solved trivially. We obtain 

x(9) = exp(rmcosh#) — 1 and x(9) = exp(rmcosh#) (43) 

for bosonic- and fermionic statistical interaction, respectively. With the help of 
these solutions we may then compute the entire scaling function. For fermionic 
statistics we obtain from (|17D 

c(r)4r m V(-ir^ l(nm) (44) 
whilst bosonic statistics yields 

6 ^ K^nrm) 

c(r) = —rm } , (45) 

7T 2 ^ n 

where Ki(x) is a modified Bessel function. We depict these functions in figure 
3(b), referring to both of them by a slight abuse of notation as A\ . One observes, 
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that the difference in the statistical interaction is effecting most severely the 
ultraviolet region. The two scaling functions converge relatively fast towards 
each other in the infrared regime. We shall encounter these feature also in other 
models. 

Using the well-known property for the asymptotic behaviour of the modified 
Bessel function lim xKAx) = 1 leads to 

x— >0 

6 1 6 

lim c(r) = -Li(—1) = - and lim c(r) = —Li(l) = 1 (46) 

for the fermionic and bosonic type equations, respectively. Here Li(x) = 



-m=l 



denotes Euler's dilogarithm [21 



We may also compute some less trivial cases. For instance taking the statis- 
tical interaction to be of "semionic" type, i.e. g = 1/2, we obtain, after solving 
(O), for the scaling function 



12r r 
c ( r ) = — T m I dOaxsvah 





- exp(— rm cosh l 



cosh 9 



= ^ m ^ n wm^^ m2n+1)rm) ■ (47) 

Once more we carry out the ultraviolet limit with the help of the asymptotics of 
the modified Bessel function 

, v 12 ~ , (2n-l)!! 12 pi* , arsinht 3 

This value of the effective central charge corresponds to the minimal model 
A^(3, 5). This is of course what we expect from section 6, since with the help 
of (|20|) it can be easily seen that we have thermodynamical equivalence between 
the Calogero- Sutherland model with coupling A = —1/2 and fermionic statistical 
interaction and the Ising model with "semionic" statistical interaction. 

For gij = 25ij we carry out a similar computation and obtain for the scaling 
function 



fir f I 1 + \/l + 4exp(— rm cosh ( 
c(r) = — m d6\n\ ¥ ^ | cosh^ . (49) 







Once again we may perform the ultraviolet limit and obtain in this case c e ff = 
2/5, which corresponds to the minimal model Ai(2,5). This is in agreement 
with the results in section 6, since the Calogero-Sutherland model with coupling 
A = 1 and fermionic statistical interaction and the Ising model with statistical 
interaction = 25ij are thermodynamically equivalent. 
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8.2 Scaling Potts and Yang-Lee Models 

Next we investigate the scaling Potts model in this context, which was previously 
studied by Zamolodchikov || with regard to conventional fermionic statistics. 
The S-matrix of the scaling Potts model equals the minimal S-matrix of y^-affine 



To da field theory [23] and reads 



sinh f § + f ) sinh f f + f ) 

Sn(0) = S 22 {9) = » and S 12 (9) = ~ J, • (50) 

smh (j - Hj smh (j - Hj 

As commented in section 4, the S-matrix does not satisfy the bootstrap at 9 = 0. 
From our point of view it seems therefore more natural to use its modification in 
the sense of (|T9| ) and employ bosonic statistics, which of course by the equivalence 
principle leads to the same TBA-equations. The two particles in the model are 
conjugate to each other, i.e. 1 = 2, and consequently the masses are the same 
mi = m 2 = m. The conjugate particle occurs as a bound state when two particles 
of the same species scatter, for instance 1+1— >2. For the TBA-equation we need 

^ = ^ = 555^1 a,ld *»W = T=£*i ' < 51) 



Then for the fermionic statistics, equation ([12]) can be solved iteratively 

In (x^(9)) = rm cosh 9 + ^ 7 d& cosh ^ z|0 ln L + _±_] 

^ (52) 

Here we assumed that the ^-symmetry of the model will be preserved such 
that x\{9) = x 2 (9). Once more for x(9) = exp(e(9)) we recover the well-known 
TBA-equation of fermionic type (first ref. in 0). It appears to be impossible 
to find analytic solutions to this equation, but it is straightforward to solve it 
numerically. Taking ln (x^^)) = rm cosh 9 one can iterate this equation as 
indicated by the superscripts. Depending on the value of mr and 9, convergence 
is achieved relatively quickly (typically n<50). The result is shown in figure 2(a) 
and appears to be in complete agreement with the calculation in 0. To make 
contact with the literature we introduced the quantity L{9) = ln (1 + x~ l {9)) . 
One observes the typical behaviour of lim mr ^ L(9) = const for some region of 9, 
which is required to derive fl34|). 

We now consider a model with fermionic statistical interaction which involves 
an S-matrix related to the scaling Potts model in the sense of fl25|) . According 
to the discussion in sections 4 and 5, this is equivalent to considering the scaling 
Potts model in which the statistical interaction is taken to be gij = 5ij — iV^. 
Since 

1 00 / 1 2 \ 

- / M<p(G) = CI] ■I A2 =[ 

3 3 



N = A--A + = -- J d9 v (9)=C~/ 2 -I A2 = [ I I ], (53) 
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(Ca 2 and Ja 2 denote here the Cartan matrix and the incidence matrix of the 
A 2 -Lie algebra), we obtain 



'-(-ID- (54) 

and the generalized TBA-equation ( ]T2|) becomes 

— oo 

(55) 

At first sight gij in ( |5*^ ) seems to be inappropriate since its off-diagonal elements 
are negative. However, due to the ^-symmetry this does not pose a problem 
(c.f. section 2), since gn + gu = and x\{6) = X2{6). In fact, this system may 



be thought of as a bosonic system with one particle species. Equation ([55]) may 
be solved easily numerically in a similar fashion as (|52]) , whereas in comparison 
with (|52D convergence is now achieved much faster (typically n < 20). The result 
is shown in figure 2(b). 

Having solved the generalized TBA-equation for x(8) we may also compute 
the entire scaling function. The result of the numerical computation is depicted in 
figure 3(b) for the two types of statistical interactions presented in this subsection. 
Notice that the conformal limit for the exotic statistics corresponds to c = 1, 
which is in agreement with fl3"T|). Hence in this limit we obtain a one particle 
bosonic system. 

Due to the fact that the scattering matrices of the scaling Potts model and 
the scaling Yang-Lee model are related as S YL (6) = S$(6)Sfe(6), the TBA- 
equations for the two models are identical. The only difference occurs due the 
fact that in the scaling Yang-Lee model there is only one instead of two particles 
present and therefore the scaling functions equal each other up to a factor 2. 

We finish the discussion of the anyonic scaling Potts model with a brief remark 
on other possible choices of exotic statistics. For example, we may take the S- 
matrix of the scaling Potts model and choose the statistical interaction to be of 
form g^ = gSij. For this choice of g^ the Z 2 -symmetry is also present, so that 



x\{9) = X2(0). This leads, in particular, to a simplification of the system fl34|) . 
In fact, it reduces to two equations of type fl38|) with i = 1, 2 and Ai = A 2 = g. 
Therefore, the effective central charge describing conformal limit of this anyonic 
version of the scaling Potts model is given by csp(g) = 2ccs(^ = g), where ccs(A) 
is the effective central charge of the Calogero-Sutherland type model discussed 
above (see figure 1). This also means that the conformal limit of the scaling Yang- 
Lee model with statistical interaction g is described by cyhig) = Ccs(A = g) (of 
course, as we discussed in section 6, this does not imply global equivalence with 
the Calogero-Sutherland model). 
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8.3 ^-minimal afRne Toda Theory 

Finally we present the fermionic computation for the easiest model, which con- 
tains at least two different mass values , i.e. y^-aftine Toda field theory. The 



minimal S-matrix reads 23 



4 ; 

v sinh(f-if) 



Hence 



*>u(*) = M*) = M*) = = (57) 

ia\ tn\ 2^2 cosh ^ 

M*) = ^(g) = - CQsh2g • (58) 

The masses are given by m x = m 3 = m/ \/2 and 777,3 = "7- Again we assume that 
the i^-symmetry is preserved, such that Xi(6) = x^{9). The TBA-equations for 
fermionic statistics read 

y y " ^2 J \ vr cosh (9-6') n cosh 2(9 - 9') 

— OO 

v v 7 Ucosh(0-0') 7rcosh2(0-0') 

—00 

where = In fl + x~ 1 (6')^ for 7 = 1,2. Once again we may solve these equa- 

tions numerically and compute the entire scaling function. The result is depicted 
in figure 3(b). The functions L^(9) exhibit the typical plateau as mr approaches 
zero. In comparison with the other models we observe that scaling functions for 
systems with the same statistical interaction have qualitatively the same shape. 
Similarly as in the previous subsection we could now also consider the anyonic 
system with = 5{j — iVy. However, in this case we would have an overall neg- 
ative contribution of p[ to the density p 2 in (f|), which as discussed in section 
2 seems inappropriate from a physical point of view. We also observed that in 
the numerical computations singularities occur. However, we may perfectly well 
choose a different type of statistics which leads to satisfactory equations. For 
instance, we may choose = gdij. In this case equations ([34]) and (|35| ) predict 
that for g = 0, 1/2, 1, 2 the conformal limit is described by c e ff taking the fol- 
lowing values: c e// (0) ss 1.16, c e// (l/2) ss 1.07, c e// (l) = 1 and c e// (2) ss 0.895. 
One may conjecture that c e ff(g) is a monotonically decreasing function. A de- 
tailed investigation of whole scaling functions for admittable choices of statistical 
interaction is left for future studies. 
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9 Conclusion 



For a multi-particle system which involves a factorizable scattering matrix de- 
scribing the dynamical interaction and a statistical interaction governed by Hal- 
dane statistics we derived the TBA-equations. These equations may be solved 
by the same means as the conventional TBA-equations of fermionic and bosonic 
type and allow the computation of the entire scaling function. 

The behaviour of the scaling functions depicted in figure 3(b) suggests the 
validity of the conjecture ||, that the series for the scaling function in the scaling 
parameter r commences with a constant and thereafter involves quadratic and 
higher powers in r, for fermionic type of statistics. For exotic statistics these 
features don't seem to be evident. 

The question of how to select a particular statistics without the prior knowl- 
edge of the nature of the particles remains to be clarified. 

It would be very interesting to generalize these kind of considerations also to 
the situation in which the dynamical scattering is described by non-diagonal S- 
matrices (e.g. affine Toda field theory with purely imaginary coupling constant), 
to the excited TBA-equations 0, Q or even more exotic statistics P3fl . 

After the completion of our manuscript Dr. Ilinski pointed out to us, that 
there exists an ongoing dispute E5j about how to achieve compatibility between 



equations (1) and (2), that is to provide a prescription for the counting of the 
states involving (2) which would lead to (1). There are even doubts whether it 
is at all possible to achieve compatibility between the two equations. What our 
manuscript concerns, this is not a crucial issue, since we could also start with a 
much weaker assumption. We need only (1) with di replaced according to (2). 

We are also grateful to Dr. Hikami for pointing out reference |2(J to us, 
in which one also finds a discussion on ideal g-on gas systems with fractional 
exclusion statistics. In particular our figure 1 is similar to the figure in there. 
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Figure 1: The effective central charge for the ultraviolet limit of the Calogero- 
Sutherland model versus its coupling constant A. 




Figure 2: Solution of the thermodynamic Bethe ansatz for the scaling Potts 
model with fermionic statistical interaction (a) and exotic statistical interaction, 
i.e. g i:j = 5ij - N ij7 figure (b). 
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Figure 3: Solution of the TBA-equation for the y4 3 -minimal affine Toda field 
theory (a). For the same value of mr the curve L 2 (9) is always below Li(6). 
Scaling function for systems involving different types of dynamical and statistical 
interaction (b). The exotic statistics for A 2 means = 8^ — iV^. 



